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Abstract 

Retarded electromagnetic potentials are derived from Maxwell's equations and 
the Lorenz condition. The difference found between these potentials and the conven- 
tional Lienard-Wiechert ones is explained by neglect, for the latter, of the motion- 
dependence of the effective charge density. The corresponding retarded fields of a 
point-like charge in both uniform and accelerated motion are compared with those 
given by the formulae of Heaviside, Feynman, Jefimenko and other authors. Con- 
trary to claims in the pedagogical literature, the fields of an accelerated charge 
given by the Feynman or Jefimenko formulae, or derived from the Lienard-Wiechert 
potentials are all different from each other as well as from those presented in the 
present paper. A mathematical error concerning partial space and time derivatives 
in the derivation of the Jefimenko equations is pointed out. 
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1 Introduction 



The present paper is the fifth in a series written recently by the present author on 
relativistic classical electrodynamics (RCED). In the first of the papers [1], all of the 
formulae of classical electromagnetism (CEM), up to relativistic corrections of 0(/? 2 ), 
relating to intercharge forces, were derived from Hamilton's Principle, assuming only 
Coulomb's inverse-square force law of electrostatics and relativistic invariance. In the 
same paper it was shown that the intercharge force, mediated by the exchange of space- 
like virtual photons, is predicted by quantum electrodynamics (QED) to be instantaneous 
in the centre-of-mass frame of the interacting charges. Recently, convincing experimental 
evidence has been obtained [2] for the non-retarded nature of 'bound' magnetic fields 
with r~ 2 dependence, (associated in QED with virtual photon exchange) in a modern 
version, probing small r values, of the Hertz experiment [3] in which the electromagnetic 
waves associated with the propagation of real photons (fields with r _1 dependence) were 
originally discovered. 

In two subsequent papers [4, 5] the predictions of the RCED formulae for intercharge 
forces derived in Ref. [1] are compared with the predictions of the CEM (Heaviside) for- 
mulae [6] for the force fields of a uniformly moving charge. Unlike the RCED formulae, 
the CEM ones correspond to a retarded interaction. If the latter are written in 'present 
time' form [7] they are found to differ from RCED formulae by terms of 0(j3 2 ). In the 
first paper [4], it is shown that consistent results for small- angle Rutherford scattering in 
different inertial frames are obtained only for RCED formulae and that stable, circular, 
Keplerian orbits of a system consisting of two equal and opposite charges are impossible 
for the case of the retarded CEM fields. The related 'Torque Paradox' of Jackson [8] is 
also resolved by use of the instantaneous RCED fields. The second paper [5] considers 
electromagnetic induction in different reference frames. Again, consistent results are ob- 
tained only in the case of RCED fields. It is demonstrated that for a particular spatial 
configuration of a simple two-charge 'magnet' the Heaviside formula for the electric field 
predicts a vanishing induction effect in the case that the 'magnet' is in motion and the 
test coil is at rest. 

In Ref. [9], the space-time transformation properties of the RCED and CEM force fields 
were studied in detail and compared with those that provide the classical description 
of the creation, propagation, and destruction of real photons. It was shown that in 
the relativistic theory longitudinal (with respect the direction of motion of the source 
charge) electric fields contain covariance-breaking terms of 0(/3 2 ). The electric Gauss 
Law and Electrodynamic (Ampere Law) Maxwell Equations are also modified by the 
addition of covariance-breaking terms of 0(/3 4 ) and 0(/3 5 ) respectively. The retarded 
fields are re-derived from the Maxwell Equations and the Lorenz condition and an error 
in the derivation of the retarded Lienard-Wiechert (LW) [10] potentials was pointed out. 
The argument leading to this conclusion — which implies that retarded fields given by the 
Heaviside formulae are erroneous for this trivial mathematical reason, as well as being 
inconsistent with QED — is recalled in Sections 2 and 3 below. 

The aim of the present paper is to present a more detailed discussion of retarded 
electromagnetic fields with a view to pointing out some of the mathematically erroneous 
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statements on this subject that have appeared in classical research literature, text books 
and modern pedagogical literature. The correct relativistic formulae for the retarded fields 
of an accelerated charge have previously been derived in Ref. [9]. These fields actually 
describe only the production and propagation of real photons whereas in text books 
and the pedagogical literature it is universally assumed that these fields describe both 
intercharge forces and radiative effects. Since the present paper is concerned only with 
the postulates and mathematical arguments used in different derivations of the retarded 
fields, the physical interpretation of the fields (in particular their relation to the quantum 
mechanical description of radiation), as discussed in Ref. [9], is not considered. 

The structure of the paper is as follows. In the following section the retarded 4-vector 
potential is derived from inhomgeneous d'Alembert equations and the Lorenz condition. 
The reason for the difference between the potential so-obtained and the pre-relativistic 
LW potentials is explained. In Section 3 Feynman's derivation of the LW potentials is 
recalled, where the 'multiple counting' committed also in the original derivations [10] is 
made particularly transparent. In Section 4 some erroneous 'relativistic' derivations of the 
LW potentials and the Heaviside formulae that are commonly presented in text books on 
classical electromagnetism are discussed. In Section 5 the retarded fields of a uniformly 
moving charge are considered and the 'present time' formulae for the retarded RCED 
fields are derived for comparison with the Heaviside formulae of CEM. In Section 6 a 
comparison is made between different formulae for the retarded fields of an accelerated 
charge that have appeared in text books and the pedagogical literature including the 
well-known ones of Feynman and Jefimenko. Section 7 contains a brief summary. 



2 Derivation of retarded electromagnetic potentials 
from inhomogeneous d'Alembert equations 



As described in Ref. [11], retarded elecromagnetic potentials may be derived from the 
Maxwell equations: 

V-E = 4ttJ (2.1) 
1 <9E 

V x B - -— = 4ttJ (2.2) 
c dt y ' 

and the Lorenz condition 

V-A+i^ = (2.3) 
c at 

where the current density J is a 4-vector: 

3(xj(t),t) = (J ; J) = (7«P*;7Wlp*) = — (2.4) 

c 

The system of source charges is assumed to be at rest in the frame S*, where the charge 
density is p*, and to move with velocity u = cf3 u relative to the frame S in which the 
potential is defined. The 4-vector velocity of the charge system in this last frame is: 

u = (c7„; cry v j3 u ) (2.5) 
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where 
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The first step of the calculation is to use the Lorenz condition (2.3) to eliminate either J 
or J from (2.1) and (2.2) to obtain the inhomogeneous d'Alembert equations: 



V 2 A - 



V 2 A 



1 d 2 Ao 
c 2 dt 2 
1 d 2 A 



-4ttJ 
-4ttJ 



(2.6) 
(2.7) 



These equations are readily solved by introducing appropriate Green functions [11]. The 
solutions give the retarded 4-vector potential: 



X n 



t) 



(2.8) 



Here x q is the position and t the time at which the potential is defined and xj(t') specifies 
the position of the volume element d 3 xj(t') at the earlier time t'. The (^-function ensures 
that the volume element lies on the backward light cone of the field point specified by 
x q , as required by causality, since the potentials give the classical description of the 
propagation, from the source to the field point, of real (on-shell) photons at speed c. 
This is a consequence of the wave-equation-like structure of the terms on the left sides of 
the d'Alembert equations. The solutions of the corresponding homogeneous d'Alembert 
equations are progressive waves with phase velocity c. 

In the special case of a single point-like source charge the current density in (2.8) is 
given by the expression: 



jQ(Sj(t'),t') ^^S(xj(t') - x Q (t')) 



(2.9) 



where XQ(t') is the position of the charge at time t'. Inserting (2.9) in (2.8), and integrating 
over xj, gives 

where 



r' = \x q -x Q (t')\, t' Q = t- 
The retarded 4-vector potential is therefore: 

( A ret. pret^ 



t- 



(2.10) 
(2.11) 



t'=t' 
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(2.12) 



A similiar result to (2.12) is obtained in the case of an extended distribution of charge 
in the case that its dimensions are much less than the separation between the average 
position of the source charge distribution, (xj), and the field point. In this case xj 
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may be replaced in the 5-function and denominator of (2.8) by (xj), so that the factor 
(r') = \x q — (xj)\ in the denominator may be taken outside the xj integral giving 



dt!_ 

V) 

dt' u 

V) 



A r f(x q ,t) = J^-J d*x J (t' J )J,(x J (t'),t')5(t' + |fg -t) 



= J ^ J d s xj(t')p%xj(t'),t')6(t' + 



l\ ,l\r/,t . \ X q { X j)\ 



t) 



where Q is the total charge of the distribution: 

Q = J p*(xj(t'),t')d 3 xj(t>) 

and 

(Q = t - 
giving the 4- vector potential: 

I 



\ 



(r f ) 



(2.13) 

(2.14) 
(2.15) 

(2.16) 



t'=(t> 



It is now of interest, in view of understanding the origin of the LW potentials, to recal- 
culate the retarded potentials after inverting the the order of the t' and xj(t') integrations 
in (2.8), so that: 

*) = / d 3 xj(t>) J dt' j0^^ 5(t' + \^~^\ _ t ) (2.17) 

Unlike in (2.10), where the insertion of the current density of a point-like charge, (2.9) 
simply specifies the value of t' in the 5-function to be t' Q , as given by Eqn(2.11), on 
integrating over xj, the argument of the 5-function in (2.17) has a more complicated 
dependence on t'\ 

5(t> - t'j) 





9 fit') 






dt' 





where t'j is the solution of the equation f(t') = and 

/(0 = t , + |fg ~ fj( * /)l ■ 

c 

It follows from (2.19), and the definition of t'j, that 



(2.18) 



(2.19) 



(2.20) 



Differentiating (2.19) gives: 



mi 

dt' 
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(2.21) 



where: 

_ x q -xj(t'j) z _ 1 dxjjt') 

rj ~ \x q -xj(t'j)y ^'c dt> {Lll) 



so that (2.17) may be written as 

AT *) = f d 3 xj(t>) f dt> J-A t* - t'j) (2.23) 

J J \x,„ — x t(t,' 1 ( 1 — r -. ■ ±. 1 



In performing the integral over proper account must now be taken of the appropriate 
current density J M to be inserted in (2.23). The limits of the t' integral are determined by 
the times at which the backward light cone of the field point coincides with the boundaries 
of the moving charge distribution. This is illustrated in Fig.l for a uniform block of charge 
DEFG, of trapeziodal shape, moving in the plane of the figure towards a distant field point, 
in this plane, far to the right. The segments AA', BB' and CC lie on the light front, 
LF, that coincides with the backward light cone of the field point. It is assumed that 
the latter is sufficiently far that LF may be approximated by a plane, with normal in the 
plane of the figure. The block of charge is moving with speed u in the plane of the figure 
at angle 9 to the direction of motion of LF. The light front starts to overlap the block of 
charge in the position AA' and ceases to do so in the position CC. The limits of the t' 
integral in (2.23) for this case then correspond to the times when the front coincides with 
AA' (lower limit) and with CC (upper limit). Inspection of Fig.l shows that, during the 
time interval between these limits, the average value of the charge density, p, is less than 
that when the distribution it at rest, p*, by the ratio: 

£ length of charge distribution 
L length of light cone overlap region 

If At' is the time during which there is overlap between LF and the block of charge, the 
geometry of Fig.l gives: 

cAt' 

L = uAt' + £ = (2.25) 

cos 

so that 

i u 

- = 1- -cos0 = l-f'j- (3 U (2.26) 
L c 

It can be seen from Fig.l that the same average charge density is obtained if the uniform 
block of charge is replaced by a point-like charge, Q, equal to the integrated charge of the 
block and placed at its centre, or if the moving uniform charge distribution is replaced 
by the fixed one MNOP with density p. For a single point-like charge the appropriate 
current density in (2.23) is then given by (2.9). (2.24) and (2.26) as: 

J«(fj(0,0 = (1 -f'j-0u)^5(xj(t') -x Q (t')) (2.27) 

Inserting (2.27) in (2.23) and performing the integrals over t' and xj recovers the result 
of Eqn(2.12). The increased overlap time of the light front resulting from the motion of 
the block of charge is exactly compensated by the reduction of the average charge density 
resulting from the same motion. The incorrect LW potentials are given by taking into 
account the time-everlap correction factor but neglecting the corresponding change in the 
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charge density. This gives, instead of (2.12), the potentials 



(A r et ;A ret ) = ( 7tt Ao(LW)^; 7u A(LW) re *) 



r'(l-r'j-p u ) 



QluPu 



t'—t' 



r'(l-r'j-p u ) 



t'=t' 



(2.28) 

where A (LW) re * and A(LW) re * are the Lienard and Wiechert potentials. This mistake 
in the original Lienard and Wiechert [10] calculations has been repeated in all text- 
book treatments of the subject of retarded potentials. Some examples may be found in 
Refs. [12, 13, 14, 15, 16]. 

Inspection of (2.23) shows that neglect of the charge density correction factor of (2.27) 
in evaluating the potentials implies that they are overestimated when the source is ap- 
proaching {f-'j ■ f3 u < 0) and underestimated when it is receding from it {f'j ■ f3 u > 0). 
On the other hand, the 1/r' dependence of the potential implies that the potentials are 

— * — * 

greater (smaller) when f'j ■ (3 U < {f'j ■ (3 U > 0). It is shown in Section 5 below that 
for the 'present time' LW fields (Eqns(4.14) and (4.15) below) the neglect of the charge 
density correction factor results in exact compensation of the 1/r' 2 dependence so that 
the magnitudes of the fields are independent of the sign of f'j ■ /3 U , as is the case for an 
instantaneous intercharge interaction. 

Note that the potentials on the right side of (2.28), derived by neglecting the density 
correction factor in (2.27), differ from the retarded LW potentials by an overall factor of 
7 n . This factor will be commented on at the end of Section 4 below where alternative 
'relativistic' derivations of the LW potentials are discussed. 



3 Feynman's derivation of the Lienard- Wiechert po- 
tentials 



The erroneous nature of the retarded potentials found when the charge density cor- 
rection factor of Eqns(2.24) and (2.26) is neglected is made particularly clear by a careful 
examination of Feynman's derivation [17] of the LW potentials for the case of parallel 
motion of the source distribution and the light front, LF, corresponding to the backward 
light cone of the field point. 

Feynman's analysis of the problem of retarded potentials is shown in Fig.2. A rect- 
angular block of charge, of uniform density, moves towards the field point, which is suf- 
ficiently far to the right that the variation of r'j may, as in deriving Eqn(2.16) above, 
be neglected in evaluating the integral that gives the potential. The light front moves 
across the charge distribution, sampling it. Each element of charge which is crossed by 
LF gives a contribution to the potential. The depth of the block of charge is I and LF 
moves over the distance L while crossing the charge distribution. The front overlaps the 
charge distribution for a time interval T. The overlap distance, L, is divided into bins 
of width w and the contribution to the potential of each bin is considered separately. In 
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Figure 1: The plane light front (LF) BB' crosses the block of charge DEFG with uniform 
charge density p while moving from the position AA' to CC'.The light front and the block 
move in the plane of the figure with speeds c and u respectively in the directions indicated. 
The average charge density sampled by the light front during its passage over the block is 
p = p*i/L The retarded potential generated by the charge of the block at a distant field 
point to the right of the figure is the same as that that would be generated by a block of 
charge in the form MNOP with the same depth as DEFG, with uniform charge density p, 
at rest, or by a moving point-like charge Q = p*V, where V is the volume of the block 
DEFG. 
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Figure 2: Feynman's method of calculating retarded potentials [17]. A uniform rectan- 
gular block of charge of length I moves to the right with speed v towards a distant field 
point. The light front, LF, in causal connection with the field point, overlaps the block 
for a distance L and a time T. In a) LF arrives at the front of the block. The position of 
the block when LF overtakes it is shown dashed, b) and d) show the positions of LF at 
times t — T/5 and t = 2T/5 respectively. The regions of the block sampled by LF in the 
time intervals < t < T/5 and T/5 < t < 2T/5 are shown by the SW-NE and NW-SE 
cross-hatched areas, respectively. The similar crossed- hatched areas in c) and e) show 
the charge volumes assigned to the potential integral, during the same time intervals, in 
Feynman's calculation. See text for discussion. 
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Fig. 2 the dimensions and velocity u are chosen so that: 



2L L 

It then follows that u = 3c/ '5. In this figure, the positions of the charge distribution 
and the front LF at times 0, T/5, 2T/5 respectively are shown. In Figs. 2b, 2d the front 
has crossed charge thicknesses of 0Aw,0.8w respectively. The region crossed during the 
time < t < T/5 is shown by SW-NE 1 diagonal cross-hatching, that crossed in T/5 < 
t < 2T/5 by NW-SE diagonal cross-hatching. Thus the average charge density in each 
bin is reduced, in comparison with the situation when the charges are at rest, by 60%. 
Integrating first over the time, as in (2.17), for each bin, then gives: 

cr j cr j 

■J bins J 

where S is the surface area of the charge distribution normal to its direction of motion 
and p is the average charge density. From the geometry of Fig. 2a, p = 2p*/5 where p* is 
the rest frame charge density. Since L = 51/2 (3.1) gives: 

A = V^! = ^ (3.2) 

cr j cr j 

where Q is the total charge in the block. Allowing for the propagation time delay of the 
light front with respect to the time of the field point (3.2) agrees with Eqn(2.16) but not 
with the LW potentials in (2.28). 

The contributions to the integral given by the first two bins, according to Feynman's 
original calculation [17] are shown by the SW-NE and NW-SE diagonal hatching in Figs. 2c 
and 2e respectively. The movement of the charge distribution is neglected, and with it the 
change in the effective charge density. Feynman's result is given by replacing p in (3.1) by 
p*, the density of the charge distribution at rest. This gives a result consistent with (2.28), 
but is evidently wrong, since charge elements are multiply counted during the passage of 
the light front. For example, a contribution to the integral is assigned proportional to 
the area of the cross-hatched region to the left of LF in Fig. 2c for t < T/5. However, 
inspection of Fig. 2b, showing the actual geometrical configuration at t — T/5, shows 
that, because of the parallel motion of the charge distribution, LF has crossed only the 
fraction of the region in Fig.2c that is both shaded and cross-hatched, not the entire cross- 
hatched region. In fact, careful inspection of Fig21-6(c) of Ref. [17] shows clearly that 
the contribution due to the passage of the light front over the first bin is overestimated. 
Only the region of the charge distribution to the left of the light front as shown in this 
figure has been sampled at this time, not the filled first bin of Fig21-6(b) of Ref. [17]. 



x The points of the compass: South- West (SW), North-East (NE), North- West (NW) and South-East 
(SE). 
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4 'Relativistic' derivations of the Lienard-Wiechert 
potentials and the electromagnetic fields of a uni- 
formly moving charge 



As well as the derivation of the LW potentials by consideration of retardation effects, 
as in the original papers of Lienard and Wiechert, text books on classical electromag- 
netism contain alternative derivations, where no retardation effects are considered, but 
instead a relativistic 'length contraction' effect is invoked. For example, in Ref. [18], the 
temporal component A of the 4-vector electromagnetic potential is obtained by Lorentz- 
transformation into the frame S, where A is defined, from the frame S* in which the 
point-like source charge Q is at rest: 

A = 7 «A* = (4.1) 

where 

I — * — * I * I — — I / a r\\ 

r=\Xg-x Q \, r =\x q -x Q \ (4.2) 

The vectors x q , xq (x*x*q) give the position of the field point and the source charge, 
respectively, in the frames S (S*). These coordinates are specified at a fixed time in the 
frame S — no retardation effects are considered. It is then assumed that the x-coordinate 
separations in the frames S and S* are related by the relativistic length contraction rela- 
tion: 

x*-x* = X q ~ XQ (4.3) 

while the y and z separations are the same in both frames. It then follows from (4.2) and 
(4.3) that 

(r * )2 = (fg ~ X Q? + ( X ~ j)i(y g - Vq) 2 + (gg ~ Z Q?\ (4 4) 

c 2 

Denoting by ip the angle between the vectors x q — xq and u, (4.4) may be written as: 

f *x2 2 [cos 2 V; + (l-g)sin 2 ^] 

r = r 



2 [l-^sin 2 ^] 



1 - 

Substituting r* from (4.5) in (4.1) than gives 



(4.5) 



A = A (LW)^ = - ® (4.6) 
r(l — p„ sin ip) 2 

This is the 'present time' (PT) formula [7] for the temporal component of the retarded LW 
potential A (LW) ret given in Eqn(2.28) above. All quantities in (4.6) are defined at the 
instant that the potential is specified. The 'present time' form of the 3-vector potential 

— * 

A is calculated, in a similar manner, to obtain 

A = A(LWr = n g^T-ri (4.7) 
r(l — sin ip) 2 
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It is interesting to note that that the ^ u factor in (4.1), manifesting the 4- vector character 
of A, is cancelled by a similar factor originating in the 'length contraction' effect of 
Eqn(4.3). 

A similar derivation of Ao(LW) PT may be found in Ref. [20] where it is noted that the 
change of variables 

X 

Xq = I Q f , V*q = Vq, Z q = Z q ( 4 -8) 

transforms the d'Alembert equation (2.6) into a Poisson equation, the solution of which is 
the Coulomb electrostatic potential Q/r*. Expressing r* in terms of (x q ,y q ,z q ), neglecting 
a multiplicative factor 7„ (which was cancelled in the derivation of Ref. [18] by the similar 
factor in the numerator of the right side of (4.1)) the potential A (LW) PT is obtained. It is 
only mentioned at the end of the calculation that the purely mathematical transformations 
of Eqns(4.8) should be interpreted as physical transformations predicted by the Lorentz 
transformation, Unlike in Ref [18] the scalar and vector potentials are treated in non- 
relativistic manner, necessitating a (tacit) neglect of a factor 7„ in order to recover the 
LW result. 

A 'relativistic' derivation of the 'present time' formulae for the electric and magnetic 
fields of a uniformly moving charge, by use of a similar 'length contraction' ansatz as in 
Refs. [18, 20] is found in Jackson's book [21] 2 , The conventional transformation laws of 
electric and magnetic fields between the frames S* and S; 

E x = El, E y = lu (E* y + p u B* z ), B z = 1u (BI + P u E;) (4.9) 

are used to transform the fields in the rest frame of the source charge: 

p * Q( x q — x q) p * Q{y q ~ Vq) p * r * r * r * n (a-\(\\ 
b * = p]3 ' = ' ^* = B * = B y = B * = ( 4 - 1U ) 

into the frame S. Performing this transformation, and using (4.3) to express the result in 
terms of S frame coordinates 3 gives 



E 



Q(x q - x Q ) 



iI{( x q - x q) 2 + (1 - $-)[(Vq ~ Vq) 2 + ( z q ~ z q) 2 }}^ 
Q cosip 

7^ 2 (1 -^sin 2 ^)i 

Q(y q - Vq) 

;/ ~ rii(x q - xq) 2 + (i - $)[(v q - vq) 2 + {*, - zq) 2 ]}* 

Q sin ip 



E 



(4.11) 



7«?" 2 (1 -Plsm ip) 



(4.12) 



B„ = P u E y = 9M^L , (4.13) 

7 2 r 2 (l -/3 2 sm ip)2 



2 A similar derivation is found in the widely-used text book on Electricity and Magnetism by Pur- 
cell [22]. 

3 Actually Jackson used a relativistic time dilatation equation equivalent to Eqn(4.3). 
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Eqns(4.11)-(4.13) may also be written in 3- vector notation as: 

E = E(H) PT = ^— - s (4-14) 

7^r 3 (l -(31 sin V) 2 

B = B(H) PT = (3 U x E (4.15) 



The label 'H' stands for 'Heaviside' who first obtained these equations [6] more than a 
decade before the advent of special relativity. They may also be obtained from the 'present 
time' potentials in (4.6) and (4.7) and the usual definitions of electric and magnetic fields 
in terms of derivatives of the 4-vector potential. 

It is easy to show that the 'length contraction' ansatz of Eqns(4.3) and (4.8) used 
to derive (4.14) and (4.15), as obtained from the retarded LW potential, but without 
invoking any retardation effect, is inconsistent with a fundamental reciprocity property of 
special relativity. This was stated in a concise way, and in a manner directly applicable 
to the problem considered here, by Pauli [23]: 

The contraction of lengths at rest in S* is equal to that of lengths at 
rest in S and observed in S* . 




Figure 3: Spatial configurations in the frames S* [a)] and S [b)] at corresponding instants 
in the two frames; for example when the origin of S*, situated at Q coincides with the 
origin of S. 



To make manifest the symmetry of the configurations in the frames S and S*, that 
is the basis of the applicability of the above reciprocity postulate in the present case, a 
test charge q, at rest, is placed at the field point in S. As shown in Fig. 3, the 'length 
at rest in S*' is the separation, r*, of the source and test charges in this frame (Fig.3a). 
Similarly the 'length at rest in S is equal to r (Fig. 3b). However, in the case of the 'length 
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contraction' ansatz of Eqns.(4.3) and (4.8), r is also the contracted value of the length r* 
as observed in S. i.e. 

r = a{u)r* (4.16) 

where ct(u) is some even function of the relative velocity u of the frames S and S*, and 
ct(u) < 1 The above reciprocity postulate states that also 

r* = a(u)r (4.17) 

Combining (4.16) and (4.17), 

r = a(u)r* = a(u) 2 r (4-18) 

It follows that if r ^ 0, a{u) 2 = 1 This is compatible with Eqns(4.3) and (4.8) only if 
u — 0, contradicting the initial hypothesis that S and S* are in relative motion. 

The length contraction ansatz of (4.3) and (4.8) is therefore incompatible with the 
above stated reciprocity property of special relativity. How this universally (until now) 
accepted length contraction effect results from a misinterpretation of the symbols in the 
space-time Lorentz transformation has been extensively discussed elsewhere [24, 25, 26]. 
In conclusion, the 'relativistic' derivation of the field equations (4.14) and (4.15) is in 
fact incompatible with special relativity and therefore fallacious. That the same result is 
obtained using the incorrect LW potentials must then be regarded, not as confirmation 
of the correctness of the formulae, but as purely fortuitous. The 'present time' formulae 
derived from the relativistically-correct retarded potentials in (2.12) are presented in the 
following section. 

An alternative 'relativistic derivation' of A(LW) PT was given by Landau and Lif- 
shitz [27]. The retardation condition (2.11) was used to write the temporal component of 
A, in the rest frame of the point-like source charge as: 

It was the noticed that the 4-vector: 

A = ( 4 -20) 
x ret ■ u 

ret ^(c(t-t' Q );x q -x Q (t')) (4.21) 

reduces to (4.19) in the rest frame of the source charge. The right side of (4.20) is 
precisely the retarded LW potential in (2.28) of a point-like charge. Although it is true 
that (4.20) gives (4.19) in the rest frame of the source charge, the same is true of the 
different 4-vector potential in (2.12). The relation (4.20) is, however, nothing more than 
a mathematical curiosity, lacking any physical motivation, whereas the potential in (2.12), 
equally consistent with (4.19), is the solution of the d'Alembert equations (derived from 
Maxwell's equations and the Lorenz condition) for a point-like charge. The physical 
meaning and method of derivation of the potential of (2.12), unlike that of (4.20), are 
therefore quite clear. 

That the retarded LW potentials and the associated fields could be derived in a 'rela- 
tivistic' calculation in which retardation effects are completely neglected, whereas in the 



where 

x 
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original derivations of Lienard, Wiechert and Heaviside, performed before the advent of 
special relativity, the (actually spurious) length contraction effect is neglected should be 
serious cause for concern. This unease, however, seemed not to be shared by authors 
of text books, and the pedagogical literature, on classical electromgnetism, throughout 
the last century. There is now ample experimental evidence for both the correctness and 
necessity of special relativity, and that retardation effects do occur in processes where 
real photons are radiated, so that the corresponding classical fields must also be retarded. 
The contradiction posed by the absence of one or the other of two essential, but different, 
physical phenomena in the two different derivations of the Heaviside formulae was clearly 
stated by Jefimenko [28], but the obvious doubt shed by this on the correctness of the 
formulae and/or the derivations, was passed over in silence. In fact, as demonstrated 
in the present paper both the original 19th century and the 20th century 'relativistic' 
derivations' are wrong. The former because the variation of the effective charge density 
of the moving charge distribution was not taken into account, the latter because the 
'length contraction' effect on which they are based, does not exist. It is proposed in the 
present paper that the correct relativistic retarded potentials of a point-like charge are 
those given above in Eqn(2.12). The corresponding electric and magnetic fields, for the 
case of a charge in uniform motion, are derived in the following section. In Section 6 the 
retarded fields of accelerated charges are considered, and compared with those derived 
from the LW potentials as well as the well-known formulae of Feynman and Jefimenko as 
well as some others that have appeared in text books and the pedagogical literature. 

The RCED 4- vector potential and current differ from those of CEM by the multiplica- 
tive factor 7 U (see Eqn(2.28) above). This leads to to a breakdown of the Gauss law for 
the electric field of a moving charge [4, 9] and covariance-breaking terms in the electro- 
dynamic Maxwell equation (Ampere's law) [9]. In text books on CEM, the validity of the 
electric field Gauss law for both static and moving source charges is justified by invoking 
the relativistic length contraction effect of Eqn(4.3) [29]. If the charge density in a moving 
frame transforms as the temporal component of a 4-vector oc 7 M , and a volume element 
transforms oc (7«) _1 due to relativistic length contraction, then the charge within the 
volume element is Lorentz invariant. Since however, as demonstrated above, the length 
contraction effect is spurious, the effective charge actually varies as 7„ or as (u/c) 2 for 
small u. This effect has been experimentally observed in the vicinity of an electrically 
neutral superconducting magnet [30]. 



Retarded electric and magnetic fields of a point- 
like charge in uniform motion 



The electric and magnetic fields corresponding to the 4-vector potential (2.12) are 
obtained by straightforward application of the definitions of electric and magnetic fields 
in terms of derivatives of the potential: 

E = -VA -i^ (5.1) 
c at 

B = V x A (5.2) 
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where, without loss of generality, it may be assumed that the electric field is confined to 
the x-y plane, 

V = i— f 

% dx q J dy q 

Unit vectors along the x- y- and z-axes are denoted as i, j and k. To perform the 
calculation, the retardation condition 



t' = t 



\x q -x Q (t')\ = t _r^_ 



(5.3) 



must be used to express the derivatives with respect to t in (5.1) in terms of t', since the 
retarded position of the source charge is a function of t', not of t. Assuming that u is 
constant, (2.12) gives: 



<9Ao 



Qlu Or' 



dx n 



Differentiating the geometrical relation: 



(r') 2 dx q 



(r') 2 =K-x Q (t')] 2 + y, 



with respect to x q gives 



, dr' 



dx n 



(x n 



Differentiating (5.3) with respect to x q , 

dt' 



OXn 



1 dr' 



c dx n 



Combining (5.6) and (5.7), rearranging, and noting that dxq(t')/dt' = cf3 u gives 

dr' 



OXr, 



t r'{l-f'-(3 u ) 



where r' = f /r'. Combining (5.4) and (5.8) 



8kr 



dx„ 



Qlu[Xq- X Q (t')] 

(l-f'-/l)(r') 3 



An analogous relation is obtained for dAo/dy q so that 

QluT 1 



-VA 



(1 - r> ■ (5 u ){r'f 

Considering now the second term on the right side of (5.1), (2.12) gives 



Idk 
c~dt 



i dk x 
c ~~dt 



c qiUq 



QluPu dr' 
{r') 2 dt' 



x qiVq 



0£ 
~dt 



(5.4) 



(5.5) 



(5.6) 



(5.7) 



(5.8) 



(5.9) 



(5.10) 



(5.11) 
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Differentiating (5.5) with respect to t'\ 

, dr' 



dt> 



= -[«,-x,(01 !fa « (0 



or 



x q>Vq 

dr' 
~dV 



df 



-cr' ■ Ax 



x q iVq 



Differentiating (5.3) with respect to t 

~di 



Cq,Vq 



! _ 1 dr[_ 
^cdi* 



x 



Combining (5.13) and (5.14) and rearranging 

~dt 



x qiVq 



1 



d^_ 
~di 



x q ->Vq 



(5.12) 
(5.13) 

(5.14) 
(5.15) 



x q ,y q 1 • (3 U 

Combining (5.1), (5.10), (5. 11), (5.13) and (5.15) gives, for the retarded RCED electric field: 



E(RCED) = 



Since A = iA, 



Similarly to (5.i 



^ 7 ? dA x 
VxA = 4 



dy q 



[r* - far* ■ fa] 
(r') 3 



j^ QluPu dr' 
(r') 2 dy q 



t'=t' 



dr' 
dy q 



So that 



B(RCED) 



ret 



V x A 



r'(l - f' ■ (3 U ) 

Qlufiu X ? 



(r'f(l-r'-p u ) 



(3 U x E(RCED) 



ret 



(5.16) 



(5.17) 



(5.18) 



(5.19) 



t'—t' 



Apart from the retarded time argument and an overall factor 1/(1 —r'-(3 u ) (the Jacobian 
of the transformation from t to t', see Eqn(5.15)) Eqns(5.16) and (5.19) are the same as 
the formulae for the instantaneous force fields of RCED [1]. 

For comparison with the Heaviside formulae (4.14) and (4.15), that may be derived 
from the LW potentials it is of interest to write E(RCED) rei and B(RCED) rei in the 
'present time' form. 

Consider a point-like charge, Q, moving with speed u along the rr-axis (Fig.4). The 
field point at which the fields are to be specified is denoted by F, the present position of 
the charge by P and the retarded position, lying on the backward light cone of F, by R. 
If N is the foot of the normal to RF passing through P, the geometry of Fig.4 gives: 

NF = r' - (3 u r' cos ip' = r'(l-f' -fa 



= \Jr 2 — (3^(r f simp') 2 = \Jr 2 — (3 2 (r simp)' 
= r (l-f3 2 u sm 2 ^=f u 



(5.20) 
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Field point 



Retarded "u Present 

position position 



Figure 4: Geometry for the calculation of the retarded potential of a charge moving with 
uniform velocity u along the x-axis in terms of the 'present time' coordinates r,ip. R is 
the position of the charge from which a light signal was emitted so as to arrive at the field 
point F at the instant shown. P is the position of the charge at this instant. The line 
segment PN is perpendicular to RF. 



Solving the quadratic equation obtained by applying the cosine rule to the triangle RFP: 

(r') 2 = r 2 + f3 2 u (r') 2 + 2(3 u rr' cos i> (5.21) 

gives the retarded separation between the source charge and field point, r' in terms of the 
'present time' parameters r and ip 

r> = r (A^t±M (5.22) 



1-0 



Also 

sin^=^= ( ]-^ ; )sin / (5.23) 

and 

,, f3 u f u + COS1p 

cosV = = - R — -r 5.24 

f' = i cos ip' + j sin vp' (5.25) 

r'-A, = p u cosij' (5.26) 

Eqns(5.20)-(5.26) may now be used to express the retarded fields in terms of 'present 
time' coordinates: 

E(RCED) PT = Q<k1 - ^MMji + cos ^ + tt l + ?A sin ^ (527) 

r 2 lu(Pu cos + f u ) 2 f u 

B(RCED) pr = ^ttt /.)»/. = * * E(RCED) PT (5,3) 
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Figure 5: E L r 2 /Q as a function of ip for various values of the source charge velocity 
f3 u = u/c. The curves of El(H)t 2 /Q are antisymmetric relative to ip = 90° and so do 
not display the expected dependence of retarded fields as seen in the El(RCED)t 2 /Q 
curves where \El(RCED)\ for ip < 90° (source charge approaching the field point) is less 
than \El(RCED)\ for ip > 90° (source charge receding from the field point). 

These expressions replace, in relativistic classical electrodynamics, the pre-relativistic 
Heaviside formulae (4.14) and (4.15). Important differences are that (5.27), unlike (4.14), 
is not radial, and in consequence does not respect Gauss' Law, and does not revert to a 
radial Coulomb field on neglecting terms of O{0£). 

The manifestly incorrect physical behaviour of the retarded electric field given by the 
Heaviside formula (4.14) is evident on comparing it with that given by (5.27). This is 
done in Figs. 5 and 6 which show curves of E£ T r 2 /Q and E PT r 2 /Q, respectively, as a 
function of ip, where E PT = £Ef T + jE^ T , for different values of f3 U) as given by (4.14) 
and (5.27). Elementary physical considerations require that if ip < 90° (i.e. the source 
charge is approaching the field point) the magnitude of the retarded field must be less 
than when if) > 90° and the charge is receding from the field point. This is because in 
the former case the source charge was further from the field point than its present-time 
position when the retarded field was emitted, and closer to it in the latter case. Because 
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5 




\|/ DEGREES 



Figure 6: Ett 2 /Q as a function of ip for various values of the source charge velocity 
j3 u = u/c. The curves of Et{H)t 2 /Q are symmetric relative to ip = 90° and so do not 
display the expected ip dependence of retarded fields as seen in the Et(RCED)t 2 /Q 
curves where \E T (RCED)\ for ip < 90° (source charge approaching the field point) is less 
than \E T (RCED)\ for ip > 90° (source charge receding from the field point). 
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the strength of the field is inversely proportional to the square of the source-field point 
separation, at the time of emission of the retarded field, the magnitude of the field must 
be greater at an angle ifj + = 90° + a than at an angle ip_ = 90° — a where a > 0. The 
fields given by (5.27) demonstrate this property, whereas E(H) PT given by (4.14) gives 
symmetric behaviour for E^: 

E T (H) PT (^ + ) = E T (H) PT (^_) (5.29) 
and antisymmetric behaviour for E^,: 

E L (R) PT (^ + ) = -E L (H)"ty_) (5.30) 

These symmetry properties are those of instantaneous [1, 4], not retarded, force fields 4 . 
As explained in Section 2 above, the antisymmetry of the El(H) pt curves in Fig. 5, 
about ip = 90° and the symmetry of the Et(H) pt curves in Fig. 6, about ip = 90° in Fig. 
6 is a consequence of the neglect of the velocity dependence of the source charge density 
in deriving the LW potentials of Eqn(2.28). 



6 Retarded electric and magnetic fields of an acceler- 
ated point-like charge: the RCED, Lienard-Wiechert. 
Feynman and Jefimenko equations 



The generalisation of the RCED formulae (5.16) and (5.19) to the case of a non- 
constant value of the source charge velocity u is straightforward. The details of the 
calculation may be found in Ref. [9]. Including the terms generated by the acceleration 
of the source charge gives: 



E(RCED) rei = 

B(RCED) re * = 
where 

It follows from (6.1) that 
A, x E(RCED) ret = 



Qifv 



i\2 



+ 



cr 



QlujPu x r') 
K 



1 llPu 



(r') 2 c/3 u r' 



t'=t' 



K = (l-f'-0 u ), P U = U 



QlujPu x r') 
K 



1 _|_ lufiuftu 



(r') 2 cr' 



QlujPu X AJ 
Kcr' 



t'=t' 



(6.1) 



(6.2) 



(6.3) 



(6.4) 



t'=t 



Q 



4 See, for example, the comparision of the 'present time' retarded LW fields with the instantaneous 
RCED fields in Figs. 2 and 3 of Ref. [4]. 
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and 



r' x E(RCED) 



ret 



QlujPu x r') 
K 



r'-Pu , llPuPu 



„A2 



+ 



cr' 



(6.5) 



t'—t 1 



The relation B(RCED) ret = p u x E(RCED) re * than holds only if $ u = 0, as in Eqn(5.19) 
above, while, in all cases, B(RCED) re * ^ f' x E(RCED) re * 

These formulae may be compared with those derived by inserting the LW potentials of 
Eqn(2.28) into the defining equations (5.1) and (5.2) of the electric and magnetic fields, 
making use of the Jacobian of (5.18) to relate derivatives with respect to t and t' . This 
calculation is given in Appendix A. It is found that: 



E(LW) ret 
B(LW) ret 



Q_ 

K 3 



r'-Pu | r'x [(r'-Pu) x p u ] 



llir') 2 

QjPu x f r 
K 3 



cr 



t'=t' 



7 2( r /)2 



+ 



- r' ■ p u ) + p u (r' ■ /3 U ) 
cr'Pu 



(6.6) 
(6.7) 



t'—t' 



The markedly different angular dependence of these fields to that of the RCED fields of 
(6.1) and (6.2) may be noticed. Eqn(6.6) gives 



P u x E(LW) 



ret 



Q_ 

K 3 



{Pu x r') 



r'-P u \ , {PuXPuW-Pu) 



+ 



cr' 



cr' 



(6.8) 



t'—t' 



P x E(LW) ret = <; 



(Pu X f') 



7 2( r /)2 



r'-P u \ | KjPu x f) 
cr' J cr' 



(6.9) 



t'—t 1 



The relations B(LW) re * = p u x E(LW) ret and B(LW) ret = f' x E(LW) re * then hold only if 
P u = 0. It is commonly incorrectly stated in text books [31, 32] that the latter relation is 
of general validity. 

The RCED retarded fields (6.1) and (6.2) are now compared with those obtained 
from formulae given by Feynman, Jefimenko and other authors. The consistency of the 
latter fields with the LW ones of (6.6) and (6.7) will also be considered. In the 'Feynman 
Lectures in Physics' formulae for the fields of an accelerated point-like charge are given, 
but not derived [33, 34]. In the notation of the present paper they are 



E(Feyn) re{ = Q 



r' 



+ 



r'd 



r' 



1 <Pf_ 
dt 2 



_(r') 2 c dt \(r')\ 
B(Feyn) rei = r'\ t , =t , x E(Feyn) ret 



+ 



t'—t' 



(6.10) 
(6.11) 



Since (see Eqn(5.3)) r' is a function of the retarded time t', not of the present time t it is 
necessary to introduce the Jacobian of Eqn(5.15) in order to evaluate the derivatives in 
Eqn(6.10). Although Feynman uses the symbol for a total time derivative rather than a 
partial one in Eqn(6.10) it is clear from the definitions of the fields in terms of potentials 
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in (5.1) and (5.2) that the time derivatives should be understood as partial ones for a 
fixed position of the field point x q as in (5.15). The straightforward but somewhat lengthy 
calculation, which is analogous that shown in the previous section, leading to Eqns(5.16) 
and (5.19) is presented in Appendix B. The following formula for the electric field is 
obtained: 



E(Feyn) 



ret 




r' + 1 [3f'(r'./3 u )-/3 tt ] 



K 



+ 



+ 



K 2 



f [3(f • AO 2 - Pl\ - 2A.(r' • A,) , [f'(f' ■ (3 U ) - fa] 



+ 



[r'(r' -Pu)-Pu] 
K 3 



(r') 



r\2 



cr' 



cr' 



t'—t' 



(6.12) 



Eqns(6.11) and (6.12) give, for the magnetic field 

Q0 U x f') 



B(Feyn) 



ret 



K 



1 2(f ' ■ (3 U ) 
Jff + K (r') 2 + K 2 



1 Uf'.p u ) 2 -(3l 



f' 



cr' 



+ 



QPu x r ' 
K 2 cr' 



(6.13) 



t'—t' 



There is no consistency between these formulae and the RCED ones of (6.1) and (6.2) or 
the ones, (6.6) and (6.7), derived from the LW potentials. One may remark the overall 
K^ 1 dependence of (6.1) and (6.2) and the overall K~ 3 dependence of (6.6) and (6.7), 
resulting from the spurious factor K~ l in the LW potentials, as compared with the more 
complicated K dependence of Feynman's equations (6.12) and (6.13). 

Contrary to what is seen above, the authors of Ref. [35] claim to demonstrate consis- 
tency of Feynman's equations (6.12) and (6.13) with fields derived from the LW potentials. 
As will be seen, no such demonstration is actually presented. The Eqn(19) of Ref. [35] is 
claimed to give the electric field derived from the LW potentials. Whether this is indeed 
the case will be discussed below. Awaiting the conclusion of this discussion, this field is 
denoted temporarily by the label LW. In the notation of the present paper it is: 



E(LW) 



t\ret 



Qr' Q d 



+ 



K(r') 2 ' cKdt' \ Kr' 
The authors of Ref. [35] then introduced the equations 



(6.14) 



t'=t' 

1 Q 



dt!_ 
~dt 



1 - 



1 dr' 




' dt 




c ~dt 


Xq 


dt' 


Xq 



K 1 c 



(6.15) 



in order to express (6.14) in terms of derivatives with respect to t (as in Feynman's 
equations (6.10) and (6.11)) instead of t' . Note that in (6.15) and the following equations 
(6.16) and (6.17) the dot denotes derivation, with respect to t, not, as elsewhere in this 
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paper, with respect to t'. Using (6.15), Eqn(6.14 ) may be written as: 



E(LW') 



A ret 




(6.16) 



The first two terms on the right side of (6.16) are the same as those of the Feynman 
formula (6.10). To obtain the latter it is therefore necessary that: 



ld 2 f' 



c dt 2 



d_ 

dt 



ur 
cr 



u 



This is not demonstrated in Ref . [35] . In fact the definition of u in terms of f is 



u 



dr* 



(6.17) 



(6.18) 



that is, in terms of a t', not of a t derivative. Substitution of (6.18) into (6.17) does not 
yield the left side of the equation. Since the last member of Eqn(21) of Ref. [35] is not 
shown to follow from the first, the purported derivation of the Feynman formulae given 
there has no basis. Indeed, as shown below, Eqn(6.14) is identical to Eqn(6.46) below as 
introduced in Ref. [37]. The calculation of Appendix A shows that Eqn(6.46) gives the 
LW electric field of Eqn(6.6). The latter is shown above (compare Eqns(6.6) and (6.12)) 
to be markedly different from the field obtained from the Feynman formula. 

A formula for the retarded potentials, similar to Eqn(2.8) above, is obtained in in 
Ref. [35] by using Green functions to solve the inhomgeneous d'Alembert equations (2.6) 
and (2.7). However, subsequently, the usual mistake is made of neglecting the motion- 
dependence of the charge density, as explained in Section 2. Thus the LW potentials 
were obtained which were claimed to give (6.14) for the retarded electric field. After 
performing the spatial integration, instead of replacing r'(t') in the argument of the 5- 
function by r'(ig), (see Eqn(2.10) above) the appropriate retarded position of the point- 
like source charge at time t, the same functional dependence on t' is assumed as before 
the spatial integration and the formula (2.18) is then used to transform the argument of 
the (^-function, leading to the spurious retardation factor l/K of the LW potentials. 

The Jefimenko formulae for the fields of an accelerated charge distribution are [36]: 



E(Jefi) ret 
B(Jefi) ret 




[p\ + 1 d[p] 



(r') 2 cr' dt 



1 ,1M 

(r') 2 cr' dt 




(6.19) 
(6.20) 



The square brackets around the charge density p and the non-relativistic current density 
J indicate that they are evaluated at the retarded time t' — t — r'/ c, as is also the spatial 
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separation r' of the current element from the field point. The volume element, d 3 xj, is also 
specified at the retarded time t'. Important differences from the RCED, LW and Feynman 
formulae are that the time derivatives act only on the charge and current densities, not 
on r', and that, as compared to the Feynman formula, only first order time derivatives 
appear. However a time derivative of f is implicit in the definition of the current J. 

In order to discuss the Jefimenko equations for the case of a point-like charge, it 
will be found convenient to explicitly impose the retardation condition by including an 
integration over the retarded time t' together with the corresponding (^-function as in 
Eqn(2.8). Indeed, much confusion about, and misinterpretation of, formulae for retarded 
fields result from not properly taking into account integrations over both space and time. 
One example is in Ref. [35] above, another will be discussed shortly. This must be done 
to correctly describe the essential physical properties of the problem under consideration 
— a spatially extended distribution of charge 5 in motion that is probed, in time, by the 
backward light cone of the field point. As will be seen, attempts to simplify formulae 
by omitting time integrals, and the associated (^-functions, as in (6.19) and (6.20), and 
in many text-book treatments of retarded potentials, often lead to erroneous results. 
Specifying explicitly the retardation condition, (6.19) and (6.20) are written as: 



E(Jefi) 



ret 




[p(xj(t'),t') , 1 d[p(xj(t'),t')} 



(r') 2 

1 d[J(xj(t'),t>)} 



+ 



cr 



dt 



c 2 r' 



dt 



S (t' + ^--t)d 3 xj(t') 



(6.21) 



B(Jefi) 



ret 



[J(xj(t'),t')} 



(r') 2 
1 d[J(xj(t'),t')\ 



cr' 



dt 



X f'5(t'+ r ^--t)d 3 xj(t') 



(6.22) 



A point-like charge Q has, in non-relativistic approximation, the following charge and 
current densities: 



p Q (x j(0,0 = QS(xj(t f ) - X Q (t')) 
J Q (xj(t'),t') = Qu(t')6(xj(t') - x Q (t')) 



(6.23) 
(6.24) 



Substituting (6.23) and (6.24) into (6.21) and (6.22) and performing the spatial integra- 
tions gives: 



E(Jefi) 



ret 



Q / dt' 



= Q 
= Q 



r' 



1 du(f) 
_(r') 2 cV dt 

du(t') 



1 dt' 



(r') 2 cV dt 



dt' 



r' 1 du(t') 

[fY ~ Kc 2 r' dt' 



5(t' - t' Q ) 



t'—t' 



(6.25) 



t'—t' 



3 In the real world consisting of an ensemble of identical point-like charged particles. 
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u(t) i du(f) 



c(r') 2 c 2 r' dt 

du(t' 



x r'S(t' - t' Q 



u(f) 1 dt' 



c(r') 2 cV dt 



dt' 



u(t') 1 du{t') 



where 



c(r') 2 Kc 2 r' dt' 

f , At' Q ) 
Z Q =t ~ — 



x r 



x r 



t'=t' r 



t'=t' 



(6.26) 



(6.27) 



For a uniformly moving charge, in contrast to the RCED, LW and Feynman equations 
the Jefimenko equations therefore predict the same (Coulombic) electric field as in the 
electrostatic case. 

In a paper on time-dependent generalisations of of the Biot and Savart and Coulomb 
laws where the Jefimenko equations were extensively discussed [37], it was claimed, con- 
trary to what is demonstrated above, that the Jefimenko, Lienard-Wiechert and Feynman 
formulae for the retarded fields are all consistent with each other. The arguments given in 
support of this conclusion are critically examined below, but first the claim of the authors 
of Ref. [37] to derive the Jefimenko equations from the defining formulae (5.1) and (5.2) 
of electric and magnetic fields and retarded potentials is considered. The assumed form 
of the potentials in Ref. [37],, in the notation and choice of units of the present paper, is: 



A 



J r' 



dr 



(6.28) 
(6.29) 



The volume element dr = d 3 xj(t') and the quantities in square brackets, as well as the 
distance r' between the volume element and the field point are all specified at the retarded 
time t' — t — r'/c. Note that unlike the solutions of the D'Alembert equations in (2.8) 
there is no integral over the retarded time in these equations and also they do not contain 
the 1/K factor of the LW potentials. Substitution of (6.28) and (6.29) into (5.1) gives 6 : 



E(GH1) 



ret 



J 



^vH + [p]v(I)+ 1 



_d[J] + [J]d_n 

c 2 r' dt c 2 dt \r' 



dr 



(6.30) 



This aleady disagrees with the corresponding equation (21) of Ref. [37], where the last 
term on the right side of (6.30) is omitted. Indeed, this term does not vanish but the last 
factor in it is: 
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(6.31) 



where (5.13) and (5.15) have been used. This result is implicit in the later Eqn(44) of 
Ref. [37], so the omission of the term in Eqn(21) of this reference in hard to understand. 



3 Thc field is labelled according to the initials of the authors, Griffiths and Hcald, of Ref. [37] 
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The retarded density [p] may be a function of the retarded time t' and the position xj(t') 
of the volume element dr, but does not depend of the position x q of the field point. 
Therefore V[p] vanishes. More formally 7 
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dx q 

at 1 



= % 



dx q 



+ . 
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dt 
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to be compared with the relation 



dx n 



i ld[pL, 
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(6.32) 



(6.33) 



given in Ref. [37]. The mathematical error leading to the incorrect equation (6.33) is a 
subtle one concerning the precise definitions of partial derivatives. Combining Eqns(5.7) 
and (5.8) gives: 
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dx r , 



1 (x q - X Q ) 

c Kr' 



so that the second line of (6.32) may be written as: 

d[p] 



dt' 



+ 



(6.34) 



(6.35) 



(6.36) 



Now it seems plausible, in view of Eqn(5.15) above to make the substitution 

1_ d_ d_ 

Kdi> ~* dt 

in (6.35) thus yielding (6.33). But all partial derivatives in (6.32) and (6.35) are evaluated 
at constant t whereas the operator relation of (6.36) is (see Eqn(5.15)) valid only at 
constant x q , and so is inapplicable in relation to derivatives with respect to x q , y q or z q . 

In fact the spurious relation (6.33) was also used by Jefimenko in the original derivation 
of Eqn(6.19) [36]. Maxwell's equations and Eqn(6.29), called the 'Vector Identity' V- 
33 [38], are used to obtain (6.19) from an integral vector identity: the 'Vector wave field 
theorem' V-31 [38]. 

The term containing V(l/r') in Eqn(6.30) is the same, up to a constant multiplicative 
factor, as one which has been previously calculated in Section 5 (Eqn(5.10)) so that: 



K{r'f 



Combining (6.30), (6.31), (6.32) and (6.37) gives: 

'W-{r' -AJ[ J] 1 d[J] 
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c 2 r' dt 
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(6.38) 



7 The ellipsis in (6.32) and subsequent equations indicates the contribution of the x- and y-components. 
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Note that the first term on the right side of (6.38) differs from the corresponding one 
in Jefimenko's formula by a factor l/K. This factor was missed in the calculation of 
Ref. [37]. In summary, the claimed-to-be-derived Jefimenko formula, Eqn(19) of Ref. [37] 
is missing a factor l/K on the first term; the second term vanishes, and the fourth term in 
(6.30) (the second in (6.38)) is omitted. Indeed, only the last term of Eqn(19) of Ref. [37] 
is correct. 



Combining (5.2) and (6.29) gives: 
B(GHl) re * = ~ c J 

Choosing [J] parallel to the x-axis 
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(6.40) 



whereas the authors of Ref. [37] state that 



v x \j ] = — y- x f' 

11 c dt 



(6.41) 



This results from a similar misuse of partial derivatives to that described above in con- 
nection with Eqn(6.33). Eqns(6.39),(6.37) and (6.40) give 
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[J] x r' 
cK(r') 2 



dr 



(6.42) 



which differs from the Jefimenko equation (6.15) by an overall factor l/K and the absence 
of the d[J]/dt term. Again the 'derivation' of the Jefimenko equation in Ref. [37] , based 
now on the incorrect formula (6.41), is fallacious. Jefimenko [36] also assumed this formula 
in order to derive the second term on the right side of (6.20). 

In Section IV of Ref. [37] it is claimed to derive the LW fields of Eqns(6.4) and (6.5) 
from the Jefimenko formulae. However this derivation starts not from the Jefimenko 
formulae (6.15) and (6.16) but instead from the different equations (given here the label 
'GH2'): 
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(6.43) 
(6.44) 



which differ from the Jefimenko equations in that the time derivatives operate not only 
on the charge and current densities but also on the reciprocal of the retarded source-field 
point separation r' and the unit vector f'. 

The authors of Ref. [37] introduce into Eqns(6.43) and (6.44) point-like non-relativistic 
charge and current densities according to Eqn(6.23) and (6.24). Since the integration over 
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t' is omitted, it is then implicit in these equations that t' = t'q, where the fixed time, t'q, 
is the solution of Eqn(6.27), corresponding to a fixed position of the source charge for 
given values of x q and t. It then follows that for point-like charges (6.43) and (6.44) are 
written as: 
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However, these formulae are not the ones obtained from (6.43) and (6.44) in Ref. [37]. 
Instead, a change of variable is introduced into the ^-functions in the first lines of (6.45) 
and (6.46): 

z(t' Q ) = x j(t' Q ) - x Q (t' Q ) (6.47) 
The Jacobian, J, relating the volume elements d 3 z and d 3 x j according to 



d 3 z = J dfx j 



3^? 



(6.4? 



is introduced. It is then stated, without derivation, that J = K where K is Jacobian 
relating dt to dt', as given by Eqns(5.15) and (6.3) above. The x-component of (6.47) is 



~Xt' Q ) = xj(t' Q ) - x Q (t' Q ) 



(6.49) 



Since xq^'q) is constant it follows from (6.49) that dz x = dxj. Similarly. dz y = dyj and 
dz z = dzj so that the Jacobian J in (6.48) is unity, not K, as claimed in Ref. [37]. 

Since the authors of Ref. [37] nevertheless did insert a factor 1/ K multiplying the 5- 
functions in the first lines of (6.45) and (6.46), before performing the spatial integrations, 
the equations obtained were not those in the last lines of (6.45) and (6.46) but instead: 
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(6.50) 
(6.51) 
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Where the subscript ' J = K' distinguishes these equations from the formally correct ones 
(6.45) and (6.46), i.e. the correctly calculated point-like charge versions of the general 
equations (6.43) and (6.44), claimed to be the Jefimenko equations but actually given, 
without derivation, in Ref. [37]. 

After writing (6.50) and (6.51) (the equivalents, in gaussian units, of the MKS Eqns(41) 
and (4.2) of Ref [37]) it is stated that they are 'essentially in the form made famous by 
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Feynman' The word 'essentially' in this context is a meaningless one. Either they are the 
same as, or they are different from, Feynman's equations on evaluating the differentials. 
This calculation, done in Appendix A below, shows that they are different. 

The introduction of the factor (1/K) inside the time derivatives in (6.50) and (6.51) is 
equivalent to replacing the retarded potentials in (6.28) and (6.29) by the LW potentials. 
It is stated in Ref [37] (again without any explicit calculation) that (6.50) and (6.51) 
yield the LW fields of Eqns(6.6) and (6.7) above. The calculation in Appendix A shows 
again that this is the case for the electric, but not for the magnetic field, since the 
latter is not given, in the general case of an accelerated source charge, by the relation 
B(LW) rei = f x E(LW) ret as claimed in Ref. [37]. 

The claim of Ref. [35] that Eqn(6.14) gives the retarded LW field will now be examined. 
Writing explicitly each term in the equation and using (5.15) to express the derivatives 
with respect to t' in terms of those with respect to t gives 



E(LW')7i* = Q 





(6.52) 



which is identical to Eqn(6.50), which is shown in Appedix A to give the LW field of (6.6). 
Thus the field given by (6.14), cannot, as claimed in Ref. [35] be the same as that given 
by the Feynman formula (compare Eqns(6.6) and (6.12) above.) 

Summarising the results obtained in this section: Ref. [35] does not demonstrate the 
consistency of the Feynman and LW formulae. The 'derivation' of the Jefimenko formulae 
from the defining equations (5.1) and (5.2) of the electric and magnetic fields and the 
retarded potentials (6.24) and (6.25) given in Ref. [37] is erroneous due to mathematical 
misinterpretation of spatial partial derivatives. The same remark applies to Jefimenko's 
original derivation [36] of these equations. The Eqns(6.39) and (6.40) given in Ref. [37] 
are not the Jefimenko equations but are obtained from them by introducing an overall 
multiplicative factor 1/K in each term and allowing the time derivatives to act on all 
factors in the terms of the equation instead of uniquely on the charge and current densities 
as in the Jefimenko formulae. Eqn(6.50) does give the LW field of (6.6), as claimed in 
Ref. [37], but the field is not 'essentially in the form of Feynman's equations (6.10) and 
(6.11). As shown in Appendix B, the electric field derived from (6.10) is markedly different 
from the LW field of (6.6). 

It was pointed out in Ref. [39] that an equation for the magnetic field identical to 
the Jefimenko equation (6.20) and a formula equivalent to the Jefimenko electric field, 
(6.19), had been given earlier in the second edition of the book 'Classical Electricity and 
Magnetism' by Panofsky and Phillips [40]. The equivalent electric field formula is 

The equivalence of Eqns(6.19) and (6.53) was shown in Section II of Ref. [39] where the 
first step was to repeat the erroneous derivation of (6.19) from the defining equation (5.1) 
of the electric field and the non-relativistic potentials (6.28) and (6.29), previously given 
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in Ref. [37] and discussed above. The term proportional to d[p]/dt is manipulated to 
obtain Eqn(6.53), As shown above, this term actually vanishes. 

7 Summary 

Retarded potentials are derived from homogeneous d'Alembert equations for electro- 
magnetic potentials and the Lorenz condition. The potentials so-obtained in Eqn(2.12) 
differ from the LW potentials of CEM. It is shown that the incorrect LW potentials result 
from neglect of the dependence of the effective density of a moving charge distribution on 
its speed. This point is made particularly clear by the careful re-examination of Feynman's 
derivation of the LW potentials presented in Section 3. In Section 4, several 'relativis- 
tic' derivations of the LW potentials or the corresponding retarded fields given in text 
books are reviewed. It is shown that they all contain misapplications of special relativity 
-in particular by invoking a spurious 'length contraction' effect. In all of the relativistic 
derivations, retardation effects are neglected, whereas in the original 19th Century deriva- 
tions of the LW potentials or the corresponding retarded fields, no relativistic effects are 
considered. There are therefore two independent, logically incompatible, and incorrect 
derivations of the LW potentials and their associated fields. In Section 5 the retarded 
RCED fields of a uniformly moving charge are derived and expressed in 'present time' 
form. Except for an overall multiplicative factor 1/(1 — r' - f3 u ) and the retarded time 
argument, they are the same as the instantaneous force fields of RCED [1]. In Section 
6 it is shown that the consistency claimed in the pedagogical literature between various 
different formulae for the fields of an accelerated charge (LW, Feynman, Jefimenko) is 
illusory. The fields all differ both from each other and from the relativistic RCED fields of 
Eqns(6.1) and (6.2). An important mathematical error concerning partial derivatives in 
the derivation of the Jefimenko formulae is pointed out. The electric field of a uniformly 
moving charge given by the Jefimenko formula is found to be, unlike CEM and RCED 
fields, Coulombic. 

The considerations of the present paper are of a primarily mathematical nature. The 
physical interpretation of retarded (radiation) and instantaneous (force) fields in RCED 
has been discussed in some detail previously [1, 9]. 
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Appendix A 



In this appendix, retarded electric and magnetic fields are derived from the LW poten- 
tials as well as from the equations (6.50) and (6.51) equivalent to those given in Ref [37] 
and claimed there to be the same as the LW fields. To derive the LW fields the potentials 
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where K = (1 — f • U ), are substituted into the defining equations (5.1) and (5.2) of 
electric and magnetic fields to give 

E(LW) rei = -VA (LW) ret - l dA(LW) ret ^ 

B(LW) re * = V x A(LW) re * (A.4) 

For simplicity, all labels, superscripts and subscripts on the fields and potentials are 
omitted in the following. 

Taking into account, by the chain rule, the contribution to the fields of each factor in 
the potentials, (A.3) and (A.4) give: 

K \r'J r' \KJ cK dt\r'J cr'dtKKJ cKr'dt 1 
B = -^xvQ-^xV^+^Vx^) (A.6) 

In these and the following equations it is understood that all spatial partial derivatives 
hold t constant and all temporal partial deivatives hold x q , the field position, constant. 
The derivatives in the successive terms on the right sides of these equations are now 
evaluated. 

The first term on the right side of (A. 5) gives: 
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where Eqn(5.8) has been used. 

Considering the second term on the right side of (A. 5), 
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Now, 
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where Eqns(5.7) and (5,8) have been used, as well as the assumption that u is parallel to 
the rr-axis. Substituting (A.9) in (A.8) and again using Eqns(5.7) and (5.8) gives: 
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The third term on the right side of (A. 5) gives 
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where (5.15) and (5.13) have been used. The fourth term on the right side of (A. 5) gives 
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where (5.13) has been used. (A.12)-(A.14) together with (5.15) then give 
QhO ( ! \ Q0u I if' • (3 U ) 2 - #) , f'-i 
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Collecting together (A. 7), (A. 10), (A. 11) and (A.15) gives, for the electric field derived 
from the LW potentials: 
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The retarded LW magnetic field given by (A. 6) is now calculated. 
Since u is assumed to be parallel to the x-axis, it follows that 
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Evaluating the first term in brackets on the right side of Eqn(A.18), 
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given by differentiating the retardation condition t' — t — r'/c as well as Eqn(5.18) have 
been used. The second tern in brackets on the right side of (A. 18) is 
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Assuming, without loss of generality, that the vector r 1 is confined to the x — y plane, 
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Combining (A.21) and (A. 22), again using (A. 20) and (5.18), gives 
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Combining A(18), (A.19) and (A.23), 
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The third term on the right side of (A. 6) is 
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where (A. 20) and (5.18) have been used. 



Collecting together (A. 17), (A.24) and (A.25), the magnetic field generated by the LW 
potentials is: 
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(A. 16) and (A.26) are the formulae (6.6) and (6.7) of the text. 



The consistency of the fields of Eqns(6.50) and (6.51) with the LW fields of (6.6) and 
(6.7) claimed in Ref. [37] is now investigated. The equations analogous to (A. 5) and (A. 6) 
given by using the chain rule to expand the derivatives in (6.50) and (6.51) are 



E = 



Qf' Qf d ( 1 \ Qf d ( 1 



Q dr' 

K(r')- cr'dt\Kj ' cK dt \r' J ~ cKr'~dt 

Q l u l) fl\ QP u d ( 1 \ Q dfi u 

cKr' dt 



cK dt \r 



or' dt \K 



5 Wu x r') Q{f3 u x f') d 

D = — — 1 - 



Kir 



A2 



cr 



dt \K 



1_\ + Q(p u x f') d 



cK dt \r 



Qf'x d(3 u 
cKr' dt 



(A.27) 
(A.28) 



34 



. Comparison with (A. 5) and (A. 6) shows that the last three terms in (A. 5) and (A. 27) 
(originating from the time derivative in (A. 3)) are the same but all other terms in (A. 27) 
and (A. 28) differ from those in (A. 5) and (A. 6). Thus to compare (A. 5) and (A. 27) 
the derivatives in the second, third and fourth terms on the right of (A. 27) must be 
evaluated, while to compare (A. 6) and (A. 28) all derivatives on the right side of (A. 28) 
must be evaluated. This is readily done using, mutatis mutandis, the formulae obtained 
above. 
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The third term is 
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Substituting (A.29)-(A.31) into (A. 22) as well as the previously obtained terms, and 
performing some algebraic simplification, gives 
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Which is the LW field of Eqn(A.16) 



Noting that the first three terms of (A. 28) differ from those of (A. 27) by the replace- 
ment f'—>j3 u x f' and using the above results for the former terms, gives, after some 
algebraic simplification, the magnetic field of (A. 28) as 
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K 3 



K 2_ R _ f3 2 f . A 



„A2 



+ 



cr' 



Q{f' x (3 U ) 



K 2 



cr 



(A.33) 



which differs from the LW field of Eqn(6.7) as well as from that given by Eqn(6.9) claimed 
in Ref. [37] to be the LW magnetic field. 
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Appendix B 



For clarity the total time derivatives in (6.6) are replaced the correponding partial 
derivatives with respect to the present time, t, for a fixed value of the field point position 



Xq . 



The second term on the right side of (6.10) contains the derivative: 



d ( f' 



dt' d 



dt' 



dt\(r') 2 J dtdt'\(r') 3 J dt 
Since the vector r 1 is confined to the x-y plane, 



1 df 1 3^ dr' 



(r') 3 dt' (r') 4 dt' 



df_ _ M x g - x q) , -fyq _ _ o 
dt' ~ % dt' +3 dt' ~ Cf u 

since dyjdt' = and cf3 u = dx Q /dt. (5.13), (5.15), (B.l) and (B.2) give: 



r' d 



c dt \(r 



A2 



1 - r' ■ P u 



3r'(r'-p u )-(3 u 
(r') 2 



(B.l) 



(B.2) 



(B.3) 



Considering now the last term on the right side of (6.10): 



d 2 f' 


dt' d 


~dt'df'~ 


dt' 


df' d 


dt 2 


dt dt' 


dt dt' _ 


dt 


at' df 



where 



d_ fdT 
dt' [~dt 



d 



d(r' • A,) 



dt> \i-r>.(3j (l-f'-(3 u ) 2 dt' 



(l-r'-f3 u ) 2 

where (A. 14) has been used. It also follows from (A. 14) that 
d 2 f' 



dt' 2 



c[f'(f' ■ (3 U ) - /3 tt ] dr' + c 



(r'Y 



dt' r' 



(B.4) 



(B.5) 



^) (r-P u )+r'(^).0 u + f'(r'J u ) 



dt 



= c 



(r') 2 



+ 



c[f'(f' ■ /i) - /t] 



(B.6) 



Combining (5.15) and (A. 14), (B.5) and (B.6) then gives 



1 d 2 f' 



(l-f'-p u ) (l-f'-(3 u y 



f'(f' ■ p u ) 2 - Pi | f'-fl, 



~A2 



cr 



+ 



r / [3(f'-/l) 2 -/3 2 ]-2/l(f'-/l) , f'(f'-p u )-p u 



„A2 



+ 



(l-r'-AO 

Inserting (B.3) and (B.7) into (6.10) yields Eqn(6.12) of the text. 



cr 



(B.7) 
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